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Abstract. A wave function of the A-component KP Hierarchy with continuous flows de¬ 
termined by an invertible matrix H is constructed from the choice of an MA-dimensional 
space of finitely-supported vector distributions. This wave function is shown to be an eigen¬ 
function for a ring of matrix differential operators in x having eigenvalues that are matrix 
functions of the spectral parameter z. If the space of distributions is invariant under left 
multiplication by H, then a matrix coefficient differential-translation operator in z is shown 
to share this eigenfunction and have an eigenvalue that is a matrix function of x. This paper 
not only generates new examples of bispectral operators, it also explores the consequences 
of non-commutativity for techniques and objects used in previous investigations. 
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1 Introduction 

The “bispectral problem” seeks to identify linear operators L acting on functions of the variable x 
and A acting on functions of the variable z such that there exists an eigenfunction z) 
satisfying the equations 

Lip = p{z)'ip and Af) = tt{x)'iI>. 

In other words, the components of the bispectral triple {L, A, 'tp{x, z)) satisfy two different eigen¬ 
value equations, but with the roles of the spacial and spectral variables exchanged. 

The search for bispectral triples was originally formulated and investigated by Duistermaat 
and Griinbaum [8] in a paper which completely resolved the question in the special case in 
which the operators were scalar differential operators with one being a Schrodinger operator. 
Since then, the bispectrality of many different sorts of operators have been considered and many 
connections to different areas of math and physics have also been discovered. (See [20] and the 
articles referenced therein.) 

The present paper will be considering a type of bispectrality in which both the operators and 
eigenvalues behave differently when acting from the left than from the right. It is necessary to 
introduce some notation and reorder the terms in the eigenvalue equations in order to properly 
describe the main results. 

Throughout the paper, M and N should be considered to be fixed (but arbitrary) natural 
numbers. In addition, LI is a fixed (but arbitrary) invertible constant N x N matrix. Let An^\ 
denote the linear functional acting from the right on functions of z by differentiating n times 
and evaluating at z = A: 

(/(z))A„,, = /(-)(A). 
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The set of linear combinations of these finitely-supported distributions with coefficients from C'^ 
will be denoted D: 


^nj,XjCj : m, {rij + 1) G N, Aj G 


Cj G C 


N 


i=i 


Finally, let W* C B) be an MA^-dimensional space of distributions. 

The goal of this paper is to produce from the selection of H and W* a bispectral triple 
(L, A, tp) such that 

• ip{x,z) = (I + is an N X N matrix function of x and z with the specified 

asymptotics in z, 

• 'ip{x, z){zH)^ is holomorphic in z and in the kernel of every element of W* (i.e., ip “satisfies 
the conditions”), 

• L'il){x,z) = 'ip{x, z)p{z) for the matrix differential operator in x acting from the left and 
some matrix function p{z), 

• and 'ip{x, z)A = 'K{x)ip{x, z) for a matrix differential-translation operator in z acting from 
the right and some matrix function 7r(x). 

In the case N = 1, this goal is already achieved constructively for any choice of distributions [23]. 
One interesting result of the present paper is that for > 1, lack of commutativity with H 
may impose an obstacle to finding such a bispectral triple in that no such triple exists for 
certain choices of W*. The three subsections below each offer some motivation for interest in 
the existence of such triples. 


1.1 New bispectral triples 

One source of interest in the present paper is simply the fact that it generates examples of bi¬ 
spectral triples that have not previously been studied. The construction outlined below produces 
many bispectral triples {L, A, ^), where L is a matrix coefficient differential operator in x, A is 
an operator in z which acts by both differentiation and translation in z, and ipi^x, z) is a matrix 
function which asymptotically approaches for a chosen invertible matrix H. 

This can be seen either as a matrix generalization of the paper [23] which considered exactly 
this sort of bispectrality in the scalar case or as a generalization of the matrix bispectrality 
in [2, 4, 32] to the a more general class of eigenfunctions and operators. 

1.2 Bispectral duality of integrable particle systems 

Among the applications found for bispectrality is its surprising role in the duality of integrable 
particle systems. Two integrable Hamiltonian systems are said to be “dual” in the sense of 
Ruijsenaars if the action-angle maps linearizing one system is simply the inverse of the action- 
angle map of the other [26]. When the Hamiltonians of the systems are quantized, the Hamil¬ 
tonian operators themselves share a common eigenfunction and form a bispectral triple [11]. 

Moreover, the duality of the classical particle systems can also be manifested through bi¬ 
spectrality in that the dynamics of the two operators in a bispectral triple under some inte¬ 
grable hierarchy can be seen to display the particle motion of the two dual systems respec¬ 
tively. This classical bispectral duality was observed first in the case of the self-duality of the 
Calogero-Moser system [21, 31]. In [23], it was conjectured that certain bispectral triples in¬ 
volving scalar operators that translate in 2 would similarly be related to the duality of the 
rational Ruijsenaars-Schneider and hyperbolic Calogero-Moser particle systems. This was later 
confirmed by Haine [19]. 
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The spin generalization of the Calogero-Moser system similarly exhibits classical bispectral 
duality [2, 32], Achieving this result essentially involved generalizing the scalar case [21, 31] 
to bispectrality for matrix coefficient differential operators. It is hoped that the construction 
presented in this paper which generalize that in [23] will similarly find application to classical 
bispectral duality in some future matrix generalization of the results in [19]. 

1.3 Non-commutative bispectral Darboux transformations 

In the original context of operators on scalar functions, the eigenvalue equations defining bi¬ 
spectrality were originally written with all operators and eigenvalues acting from the left. How¬ 
ever, research into non-commutative versions of the discrete-continuous version of the bispectral 
problem in non-commutative contexts found it necessary to have operators in different vari¬ 
ables acting from opposite sides in order to ensure that they commute [5, 14, 15, 16, 17, 18] 
(cf. [9]). Continuous-continuous bispectrality in the non-commutative context is also a subject 
of interest [6, 13, 27, 33]. As in the discrete case, it was found that the generalizing the results 
from the scalar case to the matrix case required letting the operators act from opposite sites 
[2, 4, 32]. Building on this observation, the present paper seeks to further consider the influence 
of non-commutativity on constructions and results already known for scalar bispectral triples. 

In this regard, the wave functions of the A^-component KP hierarchy are of interest since they 
asymptotically look like matrices of the form where H is an N x N matrix [3, 7]. Conse¬ 

quently, unlike the scalar case or the case H = I considered in [2, 32], the vacuum eigenfunction 
itself may not commute with the coefficients of the operators. In fact, for the purpose of more 
fully investigating the consequences of non-commutativity for bispectral Darboux transforma¬ 
tions, this paper will go beyond the standard formalism for the Al-component KP hierarchy by 
considering the case in which H is not even diagonalizable and therefore has a centralizer with 
more interesting structure. Furthermore, following the suggestion of Griinbaum [13], the present 
paper will consider the case in which both of the eigenvalues are matrix-valued. 

By generalizing the construction from [23] to the context in which the vacuum eigenfunc¬ 
tion, eigenvalues, and operator coefficients all generally fail to commute with each other, this 
investigation has identified some results that are surprisingly different than in the commutative 
case. For example, it is shown that in this context there exist rational Darboux transformations 
that do not preserve the bispectrality of the eigenfunction (see Section 6.1) and that bispectral 
triples do not always exhibit ad-nilpotency (see Remark 6.1). These will be summarized in the 
last section of the paper. 


2 Additional notation 

2.1 Distributions and matrices 

Let M, N, H and W* be as in the Introduction. The set of constant A^-component column 
vectors will be denoted by and CnxN is the set of A" x A" constant matrices. Associated to 
the selection of H one has 

C = {<5 G CnxN ■ [Q, H] = 0}, 
the centralizer of H in CtvxW- 

Let {Ji,..., 5mn} be a basis for W* C O. Unlike the selection of A, M and H which were 
indeed entirely arbitrary, two additional assumptions regarding the choice of W* will have to be 
made so that a bispectral triple may be produced from it. However, rather than making those 
assumptions here at the start, the additional assumptions will be introduced only when they 
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become necessary. This should help to clarify which results are independent of and which rely 
on the assumptions. 

Nearly all of the objects and constructions below depend on the choice of the number N, 
the matrix H and the distributions W* that have been selected and fixed above, but to avoid 
complicating the notation the dependence on these selections will not be written explicitly. (For 
instance, the matrix <1> in (3.1) could be called ^n,m,h,w* because it does depend on these 
selections, but it will simply be called <h.) 

2.2 Operators and eigenvalues 

The operators in x to be considered in this paper will all be differential operators in the variable x 
(also sometimes called ti) which are polynomials in d = having coefficients that are N x N 
matrix functions of x. The operators in ^ will be written in terms of or the translation 

operator : f(z) i— f{z + a). More generally, they will be polynomials in these having 
coefficients that are N x N matrix rational functions of z. 

Because operator coefficients and eigenvalues will be matrix-valued, the action of an operator 
will depend on whether its coefficients multiply from the right or the left. It also matters whether 
the eigenvalue acts by multiplication on the right or the left of the eigenfunction. This paper 
will adopt the convention that all operators in x that are independent of z (whether they are 
differential operators or simply functions acting by multiplication) act from the left and that 
all operators in z that are independent of x (including functions, translation operators, finitely- 
supported distributions and differential operators) act from the right. The action of an operator 
in 2 : will be denoted simply by writing the operator to the right of the function it is acting on. 
So, for instance, the function satisfies the eigenvalue equations 

The decision to have operators in x and 2 ; acting from different sides is not merely a matter of 
notation. The need for such an assumption for the differential operators in x and 2 ; respectively 
was already noted in prior work on matrix bispectrality [2, 4, 14, 15, 16, 17, 18, 32]. The present 
work extends this convention to the eigenvalues and finitely-supported distributions as well, and 
does so because the theorems fail to be true otherwise. 

Remark 2.1. Note that one needs to be cautious about applying intuition about eigenfunc¬ 
tions in a commutative setting without considering how non-commutativity may affect it. For 
example, although a non-zero multiple of an eigenfunction in the commutative setting always 
remains an eigenfunction with the same eigenvalue, here there are two other possibilities. Sup¬ 
pose Li/j^x, z) = z)p{z), so that V' is an eigenfunction for L with eigenvalue p, and that g is 

an invertible constant matrix. Then gifj may not be an eigenfunction for L if [L, 5 ] 7 ^ 0 and more 
surprisingly even though ^g is an eigenfunction for L, the corresponding eigenvalue changes 
to g~^pg. 

3 Dual construction for AT-component KP 

The purpose of this section is to produce a matrix coefficient pseudo-differential operator sat¬ 
isfying the Lax equations of the multicomponent KP hierarchy introduced by Date-Jimbo- 
Kashiwara-Miwa [7] and mostly follows the approach of SegaKWilson [29]. The proof methods 
utilized here are rather standard in the field of integrable systems. However, one of the main 
points of this paper is that some of the novel features of this situation pose unexpected obstac¬ 
les to the standard methods used to study bispectrality. So, although it is not surprising that 
the operators produced in this way satisfy these Lax equations, the proofs are presented with 
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sufficient detail to ensure that they work despite the non-diagonalizability of H and the fact 
that the distributions here are acting from the right. 


3.1 The Af-component Sato Grassmannian 

Let denote the Hilbert space of square-integrable vector-valued functions —)• 

where C C is the unit circle \z\ = 1 and is the set of complex valued row^ iV-vectors. 

Denote by e* for 0 < i < N — 1 the 1 x N matrix which has the value 1 in column i + 1 and zero 
in the others. This extends to a basis {e*: i G Z} of for which a = when i = aN + b 
for 0 < b < N — 1. The Hilbert space has the decomposition 

y^iN) =^1^ 


where is the Hilbert closure of the subspace spanned by for 0 < i and is the Hilbert 
closure of the subspace spanned by e* for i < 0. 

Definition 3.1. The Grassmannian Gr^'^^ is set of all closed subspaces V C such that the 
orthogonal projections V —)• is a compact operator and such that the orthogonal projection 

V —is Fredholm of index zero [3, 28, 29]. 

The notion of Wcomponent KP hierarchy to be considered in this paper is compatible with, 
but somewhat different from that addressed by previous authors as the following remark ex¬ 
plains. 


Remark 3.2. For the Wcomponent KP hierarchy, the construction of solutions from a point 
in the Grassmannian usually involves a collection of diagonal constant matrices Ha (1 < ct < 
N) such that powers of zHa infinitesimally generate the continuous flows and ^-dependent 
matrices {1 < /3 < N — 1) that generate discrete flows (sometimes called “Schlesinger 
transformations”) of the hierarchy [3, 7]. In the present paper, however, only the continuous 
flows generated infinitesimally by powers of z times the (not necessarily diagonal) matrix H 
selected earlier will be considered. 


3.2 A point of associated to the selection of distributions 

As usual, one associates a subspace of to the choice of W* by taking its dual in and 
multiplying on the right by the inverse of a matrix polynomial in z whose degree depends on 
the dimension of W* (cf. [21, 22, 29, 30] where the analogous procedure involved dividing by 
a scalar polynomial): 

Definition 3.3. Let W C be defined by 

W = {p(z)(zH)-^: p(z) G (p)6 = 0for 6 gW*}. 

Lemma 3.4. W G 

Proof. The image of W under the projection map onto is contained in the finite-dimen¬ 
sional subspace spanned by the basis elements for —MN < i < —1. This is sufficient to 
conclude that the projection map is compact. The map w i—)• w{zH)^ from W to 77^^ has 
Fredholm index MN because it has no kernel and the image is the common solution set oi MN- 
linearly independent conditions. The map from which first right multiplies by {zH)~^ 

^In previous papers the elements of have been written as column vectors. However, because the con¬ 

struction of bispectral operators below is most easily described in terms of matrix finitely-supported distributions 
in 2 acting from the right, they will be written here as 1 x N matrices. 
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and then projects onto has index —MN since its kernel is spanned by the basis vectors e* 
with 0<i<MA^ — 1 but the image is all of The composition of these maps is the 

projection from W to and so its index is the sum of the indices which is zero. ■ 


3.3 Af-component KP wave function 


OO 

Definition 3.5. Let V'o = exp ( ^ where t = • ■ ■) are the continuous KP time 

i=l 

variables, with the variables ti and x considered to be identical. Let </>i(t) = (1 < * < 

MN) be the C'^-valued functions obtained by applying each element of the basis of W* to V’o- 
Combine them as blocks into the N x MN matrix (f) = • • • (Pmn) and define the matrix $ 

as the MN X MN block Wronskian matrix 


ch(t) 


/ 







/ 4>i 4>2 

(t>'l </>2 


</>MV \ 

4>'mn 


AM-1) 

^MN 




{3.1) 


Assumption 3.6. Henceforth, assume that W* was chosen so that the matrix <1> in (3.1) is 
invertible for some values of x = ti (i.e., so that det<h ^ 0 ). 


The following remarks offer two different interpretations of the fact that Assumption 3.6 is 
necessary here but not for the analogous result in the scalar case [23]. 

Remark 3.7. When N = 1, the requirement that det(d>) 7 ^ 0 is equivalent to the requirement 
that {</>!,..., (Pm} is a linearly independent set of functions. Then, the independence of the basis 
of distributions would already ensure that Assumption 3.6 is satisfied. However, when A > 1 
the block Wronskian matrix can be singular even if the functions (pi are linearly independent 
as functions of x. (For example, consider the case M = \, N = 2, (pi = {1 (p 2 = (x x)~^). 


Remark 3.8. The determinant of can be interpreted as the determinant of the projection 
map from 'ipQ^{t)W to T-6^\ In other words, it is the r-function of W. (The proof of this claim 
is essentially the same as the proof of Theorem 7.5 in [22].) The r-function is non-zero when 
'ipQ^{t)W is in the “big cell” of the Grassmannian. In the case N = 1, the orbit of any point W 
in the Grassmannian under the action of ip^^ intersects the big cell [29]. In contrast, for the 
A-component KP hierarchy it is known that there are points in the Grassmannian whose orbit 
under the continuous flows never intersect the big cell [3]. 

Definition 3.9. Due to Assumption 3.6, we may define the differential operator K as 


K = d^I-[(P[^^ 


^MN 


$ 


-1 


/ I \ 

dl 


Lemma 3.10. 


(3.2) 


(o) The operator K defined in (3.2) is the unique monicNxN differential operator of order M 
such that? K(p = 0. 

^The expression Dcf) = 0 is a convenient way to write that applying the N x N differential operator D to each 
function (pi {1 < i < MN) results in the zero vector of C^. 
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(6) If L is any N x N matrix differential operator satisfying Lcjj = 0 then L = QoK for some 
differential operator Q. 

It is easy to check that Kef = 0. Alternatively, Lemma 3.10(a) follows from results of Etingof, 
Gelfand and Retakh on quasi-determinants [10]. However, Lemma 3.10(b) is apparently a new 
result. Although the lemma was originally formulated for this paper^, a self-contained proof is 
being published separately as [24]. 

Remark 3.11. The theory of quasi-determinants is not explicitly being used here but the ope¬ 
rator K defined in (3.2) could alternatively be computed as a quasi-determinant of a Wronskian 
matrix with N x N matrix entries [10]. The method of quasi-determinants was applied to the 
bispectral problem for matrix coefficient operators in [4]. So, in this sense, the use of this 
operator K here is a continuation of the approach adopted there. 

Definition 3.12. Let z) = . This function will play an important role as the 

eigenfunction for the operators in x and (given one additional assumption) z to be introduced 
below. 

Theorem 3.13. The function defined above has the following properties: 

• z) = {I + O{z~^))'if>o where I is the N x N identity matrix, 

• ipifz) G W for all t in the domain off. 

Consequently, f = fw is the N-component KP wave function of the point W G . 

Proof. Because d{fo) = f^zH, applying the monic differential operator K = -!-••• to V'o 

produces a function of the form Kfo = P{t,z)fo where P is a polynomial of degree M in z 
with leading coefficient . Then KfQH~^z~^ = {I O{z~f )fo as claimed. It remains to 
be shown that f is an element of W. By Definition 3.3 it is sufficient to note that for each 
1 < i < MN, fz^H^ = Kfo satisfies 


{Kfo)5i 


K{ifo)Si) = Kfi = 0 . 


3.4 Lax equations of the Al-component KP hierarchy 

Definition 3.14. Let K~^ denote the unique multiplicative inverse of the monic differential 
operator K in the ring of matrix-coefficient pseudo-differential operators and let C = KodoK~^ 
be the pseudo-differential operator obtained by conjugating d hy K. 

Theorem 3.15. C satisfies the Lax equations 

i-e = ncu,c] 


for each i gN, where ( ^ S cxi{t)d^. 

i=—oo 2=0 

Proof. Let f = {fo)^ for some 5 G W*. A key observation is that for each i G N: 


^Lemma 3.10(b) is used in the proofs of Theorem 3.15, Lemma 5.10 and Theorem 4.7. 
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Using the fact that </> satisfies these “dispersion relations” and the intertwining relationship 
O o K = K o we differentiate the identity K[(f) = 0 (which follows from Lemma 3.10) by ti 
to get 

0 = + K{^u) = KtM + = KtM + r o K{^) 

= KtM + (n- o m) + {c)+ o K{cp). 

Since both K and (£*)+ are ordinary differential operators (as “+” denotes the projection onto 
the subring of ordinary differential operators), (>C*)+ o AT is an ordinary differential operator 
with a right factor of K. Therefore, (p is in its kernel and the last term in the sum above is zero. 
We may therefore conclude that 

Ku{4>) + {C).oK{ct>)=^. (3.3) 

(Note that the second term in this sum need not be zero since (T*)_ is not an ordinary differential 
operator.) 

Using o K = K o we can split £* into its positive and negative parts to get 
(£*)_ oK = Kod^ - (£*)+ o K. 

Since the object on the right is just a difference of differential operators we know that (T*)_ o K 
is a differential operator. 

According to (3.3), T{cf)) = 0 where T is the ordinary differential operator 
T = Kt^ + {C)_ o K. 

Then by Lemma 3.10, there exists a differential operator Q so that T = Q o K. However, T 
has order strictly less than M since the coefficient of the order term of K is constant 
by construction and since multiplying by (>C*)_ will necessarily lower the order. This is only 
possible if Q = 0 and T is the zero operator. Hence, Kt- = —(£*)_ o K. The Lax equation 
follows because 

= Kt^o do K-^ - Kodo K-^ o o 
= -(£*)_ oKodo K-^ + Kodo K-^ o (£*)_ oKo 

= [A(n-] = [(^%,^]- ■ 

4 Operators in x having ^ as eigenfunction 

Remark 4.1. Prom this point onwards, the goal is to determine whether the wave func¬ 
tion is an eigenfunction for an operator in a: = with z-dependent eigenvalue and 

vice versa. The higher indexed time variables will only complicate the notation. So, it will 
henceforth be assumed that L = 0 for i > 2. Then, cf) and the coefficients of A can be con¬ 
sidered to be functions of only the variable x, and 'i/joix, z) = and z) (sometimes called 
the “stationary wave function”) are functions of x and z. Unlike the numbered assumptions, 
this one is made for notational simplicity only. Dependence on the KP time variables can be 
added to the objects to be discussed below so that all claims remain valid. 

Definition 4.2. A distribution <5 G D can be composed withp( 2 ;) G by defining p(z)o(5 

to be the distribution whose value when applied to f{z) is the same as that of 6 applied to the 
product f{z)p{z) for any f{z) G Ca^xAtI-^]- Associate to the choice W* of distributions the ring 
of polynomials with coefficients in C which turn elements of W* into elements of W*: 

A = {p{z) G C[z]: p{z) o 5 G W* Vh G W*]. 
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In particular, for each p ^ A and each basis element 5i there exist numbers Cj such that po 5i = 
MN 

i=l 

As in the scalar case, elements of A are stabilizers of the point in the Grassmannian: p € A 
then Wp C W. The main significance of the ring A is that there is a differential operator Lp 
of positive order satisfying Lpif; = ipp{z) for every non-constant p € A. Interestingly, unlike the 
scalar case, it will be shown that the question of which constant matrices are in A is also of 
interest in that whether ip is part of a bispectral triple is related to whether H is an element 
of A. Before those facts are established, however, the following definitions and results show 
that A contains polynomials of every sufficiently high degree. 

Definition 4.3. Let S C C denote the support of the distributions in W*. That is, A G S' if 
an only if An^\ appears with non-zero coefficient for some n in at least one 6i. For each A G S 
let mx denote the largest number n such that appears with non-zero coefficient in at least 
one element of W*. The scalar polynomial 

Po{z) = ii{z-xr^+^ 

Aes 

will be used in the next lemma, in Definition 5.8 and also in Theorem 5.13 below. 

Lemma 4.4. For any p G C[z], the product po{z)p{z) is in A. So, po{z)C[z] C A. 

Proof. Let p G C[z] and 6 G W*. Then applying the distribution popo 5 to any polynomial q 
is equal to {qpop)6. The distribution 6 will act by differentiating and evaluating at z = X for 
each A G S. However, po was constructed so that it has zeroes of sufficiently high multiplicity 
at each A to ensure that this will be equal to zero. Hence pop o (5 is the zero distribution, which 
trivially satisfies the criterion in the definition of A. M 

n 

Definition 4.5. For p G C[z] where p = Ci G C define 

i=0 


i =0 

Lemma 4.6. For anyp G C[z\, the constant coefficient differential operator \)~^{p) has ipo = e^^^ 
as an eigenfunction with eigenvalue p{z). 

Proof. If p is the polynomial with coefficients Ci as in Definition 4.5 then 


V *=0 / i =0 i=0 

= CiH-\zHy] = V-O ( x; CizA = f;op{z). 


\2 = 0 / \2 = 0 / 

Theorem 4.7. For any p G A there is an N xN ordinary differential operator Lp in x satisfying 
the intertwining relationship 


Lpo K = K o\)-^{p) 


and the eigenvalue equation Lp'if{x,z) = 'fp{z). 
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Proof. Let p G A. Observe that 

K o \>~^{p){(l)i) = K o b“^(p)((V’o)'5j) = K{{\?~^{p)tpo)Si) = K{{'tpo)p{z) o 6i) 

( MN \ MN MN 

= X] ° ^ 

i=i / i=i j=i 

But, that means that each function (f)i is in the kernel of AT o \)~^{p) and hence by Lemma 3.10 
there is a differential operator Lp such that Ko\)~^{p) = LpoK. This establishes the intertwining 
relationship. 

Applying Lp to = K'iPq{zH)~^ one finds 

Ip'll; = LpiKM^H)-^) = {Lpo = {K ob-\p))M^H)-^ 

= KipopizIizHy^ = K'tpQ{zH)~^p{z) = ijjpiz). ■ 

5 Operators in 2 : having ?/? as eigenfunction 

Definition 5.1. For any a G C let the translation operator act on functions of 2 ; according 
to the definition 

{f{z))Ta = fiz + a). 

Further let = Ylij~^a'yi{o(dzyCij where the matrices Cij and constants ji are defined by 
the formula 

exp exjp{'yix)x^Cjj (with 7 * A 1i' if * 7 ^ i')- (5-1) 

ij 

Lemma 5.2. The differential-translation operator has V’o = as an eigenfunction -with 
eigenvalue e"*.‘ 

Proof. By definition, 

gXzHjH ^ ^xzH Ta^^(adi)Cij = e^^^ ^ {axH^Cij. 

i,j id 

Then the claim follows if it can be shown that the sum in the last expression is equal to e“*. 

However, (5.1) holds not only for any scalar x but also when it is replaced by some matrix 
which commutes with the matrices H~^ and Cij which appear in it. In particular, since axH 
commutes with H~^, its commutator with the expression on the left side of equation (5.1) is 
zero. From this one can determine that its commutator with each coefficient matrix Cij on the 
right is zero as well. However, replacing x with axH yields the formula 

exp(ax) = eyip{'jiaxH){axHy Cjj 
id 

as needed. ■ 

The goal of this section is to produce operators in z which are matrix differential-translation 
operators in 2 ; having rational coefficients that share the eigenfunction if with the differential 
operators Lp in x constructed in the previous section. In order for the construction to work, an 
additional assumption is required: 
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Assumption 5.3. Henceforth, it is assumed that H € A. Equivalently, assume that for each 
1 < i < MN there exist numbers Cj such that 

MN 

H5i = ^ Cj5j. 
i=i 


5.1 The anti-isomorphism 

This section will introduce an anti-isomorphism between rings of operators in x and z respectively 
such that an operator and its image have the same action on i/jq. The use of such an anti¬ 
isomorphism as a method for studying bispectrality was pioneered in special cases in [30] and [25] 
and extended to a very general commutative context in [1]. (Additionally, three months after 
the present paper was posted and submitted, a new preprint by two of the same authors as [1] 
appeared which seeks to further generalize those results to the non-commutative context [12].) 

Definition 5.4. The two rings of operators of interest to this construction are 
W = 0e"^C[x,5] and = 0 T^C[z,9], 

CfSC aSC 

Note that operators from both rings involve polynomial coefficient matrix differential opera¬ 
tors which commute with the constant matrix H, but elements of W may also include a finite 
number of factors of the form while elements of may similarly include factors of T^ for 
a finite number of complex numbers a. 

Definition 5.5. Let b: W —)• be defined by 

( I m n \ I m n 

i=0 j=0 fe=l J i=0 j=0 k=l 

where Cijk G C are the coefficient matrices. 

Lemma 5.6. For any Lq G W, the operators Lq and \>{Lq) have the same action on '00 = e^^^: 


Ao0o — 0ob(Lo). 


Moreover, b is an anti-isomorphism. 
Proof. Since 


X] X] X] 

i=0 j=0 k 

m n 




X e 


9^(0o) = EEE 

i=0 j=0 k 


EEEg 

i=0 j=0 k 
m n 






{'lpQ){zHy = EE Eg. e<^^^{i^o){H-^dy)\zHy 

i=0 j=0 k 


i=0 j=0 k j=0 j=0 k 




(’^«)EEE(iG‘««)‘T"(.j?yc« = EEEg 


I Ofca: 


X e 




i=0 j=0 k 


i=0 j=0 k 


it follows that Lq and b(Lo) have the same action on 0o- 
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The inverse map b ^ is found by simply moving the factor of * to the other side of 
equation (5.2), which confirms that b is a bijection. If Ki, K 2 G W we have 

{K2 o = K 2 {Ki{i^o)) = i^2(V’ob(i^i)) = (A:2(V'o))b(i^i) 

= (V’ob(i^2))b(i^i) = V'o(b(i^2) o b(Ai)). 

However, we also know that {K2 o = il^ob{K2 o Ki). Then, b(Ar 2 o Ki) — \>{K2) o \>{Ki) 

has the V’o in its kernel. The only operator in having ipo in its kernel is the zero operator. 
Therefore, 

b(A 2 o Ki) = b{K2) o b(i^i) (5.3) 

and the map is an anti-isomorphism. ■ 

Remark 5.7. The reader may be surprised to see “anti-isomorphism” being used to describe 
a map satisfying equation (5.3). Generally, a map having this property is called an isomorphism. 
There are two reasons this terminology is being used here. First, that is the terminology that was 
used for the analogous map in previous papers in the scalar setting [1, 23, 25]. More importantly, 
the fact that the order of operators here is preserved is merely a consequence of the fact that 
the operators in z act from the right while the operators in x act from the left. Note that when 
K 2 o Ki acts on a function it is Ki that acts first while when b(A' 2 ) o b(iFi) acts it is b(Ar 2 ) that 
acts first and it is in this sense that it is an anti-isomorphism. Nevertheless, it is interesting 
to note that in this context the map does preserve the order of operators in a product. It 
may be that had the operators in different variables been written as acting from opposite sides 
even in the commutative case from the beginning, the map would have been described as an 
isomorphism instead. 

5.2 Eigenvalue equations for operators in 2 

Definition 5.8. Let Lq be the constant coefficient matrix differential operator 
Lo = b~^{po{z)I). 

Lemma 5.9. LqcJ) = 0 and so Lq = Q o K for some Q. 

Proof. Because differentiation in x and multiplication on the left commute with the application 
of the distributions <5* we have 

Lo(l)i = Lo{'il>o)5i = {Lo^po)6i. 

By Lemma 5.6, Lq^jq = 'iPqPq{z). The polynomial po was chosen so that it has a zero of high 
enough multiplicity at each point in the support of 6i to guarantee that 


{^poPo{z))di = 0 . 


It then follows from Lemma 3.10 that Lq = Q o K. ■ 

In the remainder of this construction, Q will denote the differential operator such that Lq = 
QoK. 

Lemma 5.10. Given Assumption 5.3, the operators Lq, K and Q commute with H. 
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Proof. Applying the operator H ^ KH to (pi we see that 

( MN \ 

( MN \ 

Y^cjKiPj)] = 0 . 

However, by Lemma 3.10, K is the unique monic operator of order M having all the basis 
functions (pi in its kernel so H~^KH = K. 

We also know that Lq commutes with H because po is a scalar multiple of the identity and 
which turns pq into Lq only introduces additional powers of H. Then 


HQoK = HLq = LqH = Qo KH = Qo HK = QH o K. 

Multiplying this by K~^ on the right yields HQ = QH. ■ 

Note that K and Q are probably not in W, as they may be rational functions in x and a finite 
number of functions of the form e“*. However, it is possible to clear their denominators either 
by multiplying by a function on the left or by composing with a function on the right, which 
motivates the following definition: 

Definition 5.11. Let be the set of x-dependent N x N matrix functions defined as follows: 

>1^ = < vr G : vr = ttqttx, 'Kk{x)K G W, Q o 'Kk{x) G W 

I aeC 

In other words, it is the set of zero order elements of W which factor as a product such the right 
factor times K and Q composed with the left factor are both elements of W. 

Lemma 5.12. jiP is non-empty and contains matrix functions that are non-constant in x. 

Proof. Note that K and Q are operators that are rational in x and a finite number of terms of 
the form e^. If we let ttk be the least common multiple of the denominators of the coefficients 
in K then is in W (since we have simply cleared the denominator by multiplication). 

Similarly, if we let ttq be a high enough power of the least common multiple of the denominators 
of Q then Q o ttq G W. Then, tt = ttqttk is by construction an element of A^. Since, 'Kq/ttk 
is also an element of for any order zero element of W, A^ contains non-constant matrix 
functions. ■ 



The main result is the construction of an operator A in z with eigenvalue vr when applied to 
the wave function ip for every vr G .A^: 

Theorem 5.13. Let tt = hq-kk £ where K = tikK G W and Q = Q o ttq £ W. Define 
A := (zH)^ o [pq{z))~^ o \>{Q) o b{K) o {zH)~^, then 

pjA = 7r(x)V’. 

Proof. Write Lq as Lq = Qo{7TK{x)7rQ{x))~^oK. Applying this operator to TpQ and multiplying 
each side by Pq^{z) gives 


{Qo {7rK{x)TTQ{x)) ^ O K)fiQPQ^{z) =1pQ. 
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Moving K to the other side using the anti-isomorphism and applying b(Q) to both sides this 
becomes 

(Q O {■KK{x)'KQ{x))~^)'iiJo{\){K) O Pq^{z) o b(Q)) = ^q\){Q). 

Moving the last expression to the other side of the equality and moving \){Q) in it to the other 
side, we finally get 

Q((7rii'vrQ)“ Vo(b(K) o p~^ o b(Q)) - V’o) = 0 . 

Note that Q is a differential operator in x with a non-singular leading coefficient (since it is 
a factor of the monic differential operator Lq). Hence, its kernel is finite-dimensional. The only 
way the expression to which it is applied, a polynomial in z multiplied by , could be in its 
kernel for all z is if it is equal to zero. From this, we conclude that 

b(A) o b(0) = b(7ri^7rQ). 

Using this one finds that the action of A on -i/; is 

(V^)A = = {7r],\x)KMzH)-^)A 

= {^T];^{x){^po)\’K{zH)~^)A = (7r^^(x)V'o)bK o {zH)~^ o A 
= {-k^^{x)'4)o)\)K o (po(^))“^ o b(Q) o \){K) o {zH)~^ 

= (7r^^(x)V'o)Ao o b(A) o {zH)~^ = {7r]^^{x)7rK{x)7rQ{x)'ilJo) o \){K) o {zH)~^ 

= {'Kq{x)'iPo) o \>{K) o {zH)~^ = (7rQ(x)7ri^(x)7r^^(x)V’o) o o {zH)~^ 

= t:{x){-kJ}{x)'iPo)\){K) o {zH)~^ = t:{x){ttJ>{x)KiPo){zH)~^ 

= 7r{x) {K'ij;Q){zH)~^ = 7r{x)'ip. ■ 

5.3 Ad-nilpotency 

We now have eigenvalue equations Lpif) = 'ipp{z) and ijjA = Tt(x)'ijj. In the commutative case, 
Duistermaat and Griinbaum noticed that these operators and eigenvalues satisfied an interesting 
relationship that goes by the name of “ad-nilpotency” [8]. As usual, for operators R and P we 
recursively define ad^P by 

ad\P = Ro P - Po R and ad^P = Ro (ad"“^P) - (ad"“^P) o P, n > 2. 

When P is a differential operator of order greater than one, then one generally expects the order 
of ad^P to get large as n goes to infinity, but they found that when P is a scalar bispectral 
differential operator and P is the eigenvalue of a corresponding differential operator in the 
spectral variable, then surprisingly ad^P is the zero operator for large enough n. In fact, in 
that scalar Schrodinger operator case they considered, Duistermaat and Griinbaum found that 
ad-nilpotency was both a necessary and sufficient condition for bispectrality [8]. 

As it turns out, in the this new context even the easier of these two statements fails to hold. 
In particular, the ad-nilpotency in the commutative case is a consequence of the fact that the 
leading coefficients of P o P and P o P are equal, but this is not generally true for differential 
operators with matrix coefficients. Gonsequently, it is not the case that ad-nilpotency holds for 
all of the bispectral operators produced by the procedure described above. The following results 
(and an example in Remark 6.1) show how and to what extent the old result generalizes. 

Theorem 5.14. Given Lp, A, n{x) and p{z) as above, the results of applying ad^Lp and ad^p 
on Tp are equal for every n G N. 
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Proof. The case n = 1 is obtained by subtracting the equations 
'K{x)Lpijj = tt{x)'4jp{z) = ijjAp{z) 
from the equations with the orders of the operators reversed 
LpTf{x)'il> = LpijjA = ijjp{z)A. 

If one assumes the claim is true for n = k then the case re = A: + 1 is proved similarly and the 
general case follows by induction. ■ 

So, the equivalence of the actions of the two operators generated by iterating “ad” remains 
true regardless of non-commutativity. One cannot conclude from this alone that either of them 
is zero without making further assumptions, but if ad^L = 0 (which would be the case, for 
instance, if one could be certain that tt would commute with the coefficients) then because V’ is 
not in the kernel of any non-zero operator in 2 : alone, the one could conclude the same is true 
for the corresponding operator written in terms of A and p: 

Corollary 5.15. If ad^L = 0 then ad\p = 0. 

Remark 5.16. This may be the first time that ad-nilpotency has been considered for translation 
operators as well as for matrix operators. It is therefore relevant to note that Corollary 5.15 is 
valid even when the operator A involves shift operators of the form Tq, (see Remark 6.1). 


6 Examples 

6.1 A wave function that is not part of a bispectral triple 

Consider the case W* = span{5i, 52 }) 


H = 


I 0 


= 


_ 


,0 -V’ V^o,o 

Then Assumption 3.6 is met and 


and 


52 = 


Ai,, 


I,.. 

x'^z xz 


\ 


The matrix polynomial p{z) = z^I satisfies po 6i = 0 for i = 1 and i = 2 and so p G A and 
as predicted by Theorem 4.7, the corresponding operator 


Ln — 




x^ 


4 


satisfies Lpip = ipp. 

Note that this is a rational Darboux transformation of the first operator and eigenfunction 
from the bispectral triple (9^, <9^, e®^^). Since in the scalar case it has been found that rational 
Darboux transformations preserve bispectrality, one might expect Lp and ip to be part of a bi¬ 
spectral triple. However, H6i 0 W* and so Assumption 5.3 is not satisfied. Thus, Theorem 5.13 
does not guarantee the existence of a differential-translation operator A in 2 ; having ip as an 
eigenfunction. In fact, in this simple case we can see that no such operator exists. 
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Briefly, the argument is as follows. Consider a differential-translation operator A acting 
on ^ from the right and suppose there is a function 7r(x) such that 'ipA = 7 r(x)'ijj. By noting 
the coefficients of and in the top right entry of each side of this equality, we see 
immediately that A 12 = 7ri2 = 0 (i.e., A and vr are both lower triangular.). Then,we have the 
scalar eigenvalue equation 'ij^nAn = vrn'i/’ii and also ?/> 2 iAii = 7r2i'0ii+7r22V'2i- Combining these 
with 'i/iii +x'ip 2 i = one concludes that e^^An can be written in the form (/(x) + g{x) jz)e^^. 
It follows that All as an operator in 2 ; has only coefficients that are constant or are a constant 
multiplied by I/ 2 ;. In fact, all of the operators in 2 ; having "01 1 as eigenfunction are known; they 
are the operators that intertwine by 5 — ^ with a constant coefficient operator. The only ones 
that meet both the criteria of the previous two sentences are the order zero operators. In other 
words. All would have to be a number. Then the equation for the action of An on ^21 tehs us 
that 7r2i = 0 and 7^22 = is also a number. Since the eigenvalue of this operator A does not 
depend on x, the operator does not form a bispectral triple with L and tp. 


6.2 A rational example with non-diagonalizable H 

Consider H = I + U, where 


U = 


0 1 
0 0 


so that 


= 


0 


The distributions chosen are W* = span{(5i,..., < 54 } with 


= ' 0 


<52 = 


0 

A2,1 


<5,2 = 


Ao,i 

0 


and 


54 = 


0 

Ao,i 


Assumption 3.6 is met and the unique monic operator of order 2 satisfying K((j)i) = 0 where 
(l>i = (V'o)<5i for 1 < z < 4 is 


(-2x - 1 


K = d^ I + 


X 

0 


X 


d + 


/x + 1 + x' 

2 


V 


X 

0 


x^ 

X -I- 1 


Now define 


/ xz‘^ — 2xz — z + X + 1 


'll; = K{'4;o){zH) = 


\ 


xz^ 


0 


z - 1 


xz^ 


^xzH 


xz^ — 2xz — z -|- X + 1 


xz^ 


Assumption 5.3 is met because H6i = 6i, H6i+i = 5i + 5i+i for i = 1,3. Consequently, it will 
be possible to produce operators in z sharing the eigenfunction ip. The matrix polynomial 

p{z) = ^ Jjs) = (^ - 1)'^ + (^ - ^ fU 

is in A because po 6i = 25j+2 and p o 5j+2 = 0 for i = 1, 2. So, 

,_i, ... (d‘^ — 2d+l — 59^ + 5(9 — 1 

*’ W"» = ( 0 #-28 + 1 


satisfies the intertwining relationship 
K o b“^(p(z)) = LpO K 
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with 


Lp — 


0 1 
0 0 


3^ + 


1 -5 
0 1 




-2 5-^ 




d + 


. 7 3 


V 


0 


1 - 


x^ 


0 -2 

This operator has eigenfunction ■0 with eigenvalue p: 

Lpi) = Lpo A'V’o = Ko LqiPq = K'ipop{z) = ippiz). 

(As mentioned in Remark 2.1, multiplying ■0 by a matrix on the right has the effect of conjugating 
the eigenvalue. It is therefore worth noting that the eigenvalue p here is not only non-diagonal 
but in fact non-diagonalizable.) 

To produce an operator in z sharing -0 as eigenfunction, we consider the polynomial poiz) = 
{z — 1)^ (which has this form because A = 1 is the only point in the support of the distributions 
in W* and because the highest derivative they take there is m\ = 2). Then 

Lq = b-^{po{z)I) = {d^ - 35^ + 3a - l)/+ (-35^ + - 3d)U 

is the operator satisfying LqV'o = V'oPo(-s) and Lq factors as Lq = Q o K with 


/a + i^ -3a+ 


Q = 


\ 


0 


a + 


X 

1 — X 


The next step is to choose two functions from C[x], so that K = ttkLC and Q = Q o ttq are 
in C[x,a]. It turns out that the selection 


ttk = 


X X 

0 X 


= '^Q 


works and we get that 


A = {zHY o {po{z)) ^ o b{Q) o b{K) o (zH) ^ = dfl -|- d: 


-2 



+ dz 


( 4 

2 {z'^ 8z -h 6) ^ 


/ 

-2z‘^ + Az + 2 

4z^ - 6z^ -8z- 20\ 

z — z'^ 

{z — 1)2^2 

+ 


{z — 1 ) 22:2 

(z — lYz^ 


4 

V 

0 

-2z^ + 4:Z + 2 

1 

to 


{z — l)^z^ / 


does indeed satisfy 'i/’A = 'Kq{x)t:k{x)Y- 


6.3 Computing for non-diagonal H 

Suppose H = \I + U where U is still the upper-triangular matrix from the previous example. 
Since 

exp(xR-^) = e^/^I - 

the corresponding operator formed by replacing with and other occurrences of x 

with adz would be 

ct a/X^ ^2 (y/X^^zU' 

This operator in z that combines differentiation and translation has the property that (e*^^)T^ 

_ ^ax^xzH 
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6.4 An exponential example with H = I 

Finally, consider H = I and 


w* / ( ^^0,0 + Ao,i\ ( 

W =span|(^ 


0 

-Ao,o + Ao,i 


Then 


K = 


id--^ 

1 + e^ 


and 


d + 


'ip{x,z) = 


V - 1 “ ' 1 - / 

{e^{z-l) + z) 

(1 + e^)z 


^zx-\-x 


\ (—1 + z 


0 ^ 

^xz (^gX^^ _ 2) — 2;) 

(-l + e^)z / 


Because of the support of the distributions, po{z) = z"^ — z and we may choose p G A to be 
a multiple of that by a constant matrix 


p{z) = 


0 - z 


z^ — z 


0 


The constant operator b ^(p) = p(9), satisfies the intertwining relationship K o p{d) = LpO K 
with 


Lp — 


/ (-1+ - 2e2^) 

(-1 + e‘^^f 

(-3 + 2e^) 


2e 


2x 


\ 


V (-l + e2^)^ 

/ 


+ 


-1 + e2^ 


-1 + 2e"^ + 2e2"^ - 2e3^ - 


V 


(-1 + e2^)^ 


(-1 + e^)^ (1 + e^) 


(-l + e*)^(l + e^) / 

-l-e"^ + 3e2*-e3^\ 
(-1 + e^)^ (1 + e^) 

(-l + e^)(l + e^) J 


d + 


1 0 / 


Hence Lpip = ippiz). (Since p{z) and V' do not commute, it is necessary to write the eigenvalue 
on the right rather than the left: Lip 7 ^ p{z)p^.) 

Now, to compute an operator in 2 : having the same function ip as an eigenfunction, we factor 
Po{d)I as Q o K to obtain 


Q — dl + g4 


0 

1 


l-e2^ -1+e^ 

Letting ttk and ttq be 


vrx = 


- 1 
0 


0 

1 -e2 


2x I ^ TTq — 


we get 


- 1 


g2x _ 


K = ttkK = 


-e^(-l + e^) 


—e 


e* + e 


2x 


+ 


(-l + e^)(l + e^) 0 

0 1 - e2^ 


dew 
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and 


Q = Q o-kq = 


1-e^ + 2e 


2x 


—e 


0 

l + e^ + 2e 


2x 


+ 


-1 + e^^ 0 

0 -1-he^^ 


dew. 


Replacing by and 5 by z we get the corresponding translation operators 


b(K) = 


0 


+ Ti 


1 0 

-1 1 


+ T2 


z - 1 0 

0 1-z 


and 


HQ) = 


l-z 


0 


+ Ti 


0 l-z 
As predicted by Theorem 5.13, 


A=^ob(Q)ob(K)o- = 
Po(z) z 


-1 0 

-1 1 


+ T2 


z + 2 0 

0 z + 2 


^ z^ + + 6z 

z{z + 2) (z -|- 3) 


V 


0 


—z^ — 5z^ — 6z 


0 

+ Ti I 2z^ + 6z-l-4 


+ T2 


+ T3 


+ T 4 


z(z + l){z + 2) 

/ 4^ + 12 

z(z -h 2)(z -h 3) 
-2z^ -4z-2 


z(z -h 2)(z -h 3) / 
l'-2z^ - 10^2 - I2z 

+ 2) (z 4- 3) 

-2z-4 


0 


0 

4z + 12 


\ z{z + 1) (z -|- 2) 


2z^ + lOz^ + 12z 
z(z 4- 2)(z 4- 3) ) 


\z{z 4-l){z 4-2) z{z 4-2)[z 4-3) / 

(z^ + hz^ + 2z-3 

z{z 4- 2)[z 4- 3) 

-z^ - bz"^ - 2z 4- • 


0 


V 


0 


z(z 4- ‘2){z + 3) / 


satisfies ip A = itip (vr = ttqttk). 


Remark 6.1. This is a good example for demonstrating ad-nilpotency and how it has changed 
in this non-commutative context. Using L, A, p{z) and ttk, ttq as above, it is indeed true that 


adlLiP = iPadlp 


(with the operator in z acting from the right as usual). However, contrary to our expectation 
from the commutative case in which the order of the operator on the left decreases to zero when 
n gets large, ad^lL is an operator of order 2 for every n. This happens because the action of 
iterating ad* on the leading coefficient of L itself is not nilpotent. On the other hand, if instead 
of t:k and ttq we had chosen 

f „‘ix 1 r ~ _ i „‘2x 1 \ 2 7- 

~ Ij-') ^ — \A ^ 

(intentionally selected so as to commute with the coefficients of any operator) then Theorem 5.13 
would have produced a different operator A* satisfying Tr*ip = ip A*. In this case, because the 
order is lowered at each iteration, ad|*L = 0. This is not particularly surprising or special as 
the same could be said if tt* were replaced by any function of the form f{x)I. However, because 
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of the correspondence in Theorem 5.14, we can conclude from this that ad^*p(z) = 0 also. That 
is more interesting because in general repeatedly taking the commutator of the operator A* 
with some function will not produce the zero operator, even if all coefficients are assumed to 
commute. Nevertheless, with these specific choices everything cancels out leaving exactly zero 
when adA* is applied three times to the function p{z) that is the eigenvalue above. 

7 Conclusions and remarks 

Given a choice of an invertible N x N matrix H and MA^-dimensional space W* of vector¬ 
valued finitely supported distributions, this paper sought to produce a bispectral triple (L, A, ip) 
where L is a differential operator acting from the left, A is a differential-translation operator 
acting from the right and ip is a common eigenfunction that is asymptotically of the form 
satisfying the “conditions” generated by the distributions. In the scalar case, this was achieved 
in [23] for any choice of W*. In the matrix generalization above, however, the construction 
only works given Assumptions 3.6 and 5.3. The bispectral triples produced given those two 
assumptions include many new examples both in the form of the eigenfunction (asymptotically 
equal to and the fact that the matrix-coefficient operator A may involve translation in z 

as well as differentiation in z. More importantly, this investigation yielded some observations 
that may be useful in future studies of bispectrality in a non-commutative context. 

This paper sought to develop a general construction of bispectral triples with matrix-valued 
eigenvalues and also to understand what obstructions there might be to generalizing the con¬ 
struction from [23] to the matrix case. It is interesting to note that these seemingly separate 
goals both turn out to depend on the non-commutativity of the ring A of functions that stabilize 
the point in the Grassmannian. Since A is also the ring of eigenvalues for the operators in x, 
it is not a surprise that by letting its elements be matrix-valued gives us a non-commutative 
ring {Lp-. p G A} of operators sharing the eigenfunction ip. It was not clear at first that the 
ring A would also be the source of the obstruction to producing bispectral triples. However, 
the construction of operators in z sharing ip as eigenfunction uses the assumption that H is an 
element of A. (Specifically, this is used in the proof of Lemma 5.10.) It is interesting to note 
that this also depends on the fact that we considered a matrix-valued stabilizer ring. 

Section 5.3 explored the extent to which the property of ad-nilpotency, which has been 
a feature of papers on the bispectral problem since it was first noted in [8], continues to apply 
in the case of matrix coefficient operators. It is still the case that the operator formed by 
iterating the adjoint of one of the operators in a bispectral triple on the eigenvalue of the other 
operator has the same action on the eigenfunction as the operator formed by iterating the adjoint 
action of its eigenvalue on the other operator. However, unlike the scalar case, only if additional 
assumptions about the coefficients of the operators are met will either of those be zero for a large 
enough iterations. 

It was previously observed [2, 4, 32] (see also [14, 15, 16, 17, 18]) that requiring the operators 
in X and z to act on the eigenfunctions from opposite side resulted in a form of bispectrality whose 
structure and applications more closely resembled that in the scalar case. Because generalizing 
the scalar results of [23] necessitated also requiring that the eigenvalues and finitely-supported 
distributions act from the same side as the operators acting in the same variables, this previous 
observation can now be extended to those other objects as well. 

This is the first time that the bispectral anti-isomorphism was used to construct bispectral 
operators in a context involving operators with matrix coefficients, see section 5.1 (cf. [12]). 
Some modifications were necessary since the rather general construction in [1] assumed that 
there were no zero-divisors so that formal inverses could be introduced. In addition, the use of 
the method here depended on the convention of considering operators in x and z to be acting 
from opposite sides and required that the coefficients of the operators on which it acted were 
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taken from the centralizer of H. The most interesting difference may have been that because 
the operators in z are acting from the right rather than the left, the map actually preserves the 
order of a product. 

Unlike the scalar case, not every choice of distributions W* corresponds to a bispectral triple. 
For instance, if Assumption 3.6 fails to be met then there simply is no wave function ijj satisfying 
the conditions. More interestingly, if Assumption 3.6 is met but Assumption 5.3 is not then there 
is a wave function that is an eigenfunction for a ring of differential operators in x, but Theo¬ 
rem 5.13 does not produce a corresponding operator in z. In fact, as Section 6.1 demonstrates, 
in at least some cases there actually is no bispectral triple of the form considered above which in¬ 
cludes that wave function. This is very different than the scalar case. Unfortunately, this paper 
does not entirely answer the question of which choices of matrix H and distributions W* produce 
a wave function V’ that is part of a bispectral triple of the type considered here. In particular, 
this paper does not show or claim that there cannot be a differential-difference operator A in z 
having ijj as an eigenfunction with x-dependent eigenvalue when H ^ A. 

Arguably, some of the non-commutativity involved in the construction above was “artificially 
inserted” in the form of the choice of the matrix H. If one chooses to consider only the case 
H = I, then Assumption 5.3 is automatically met and Theorems 4.7 and 5.13 produce a bi¬ 
spectral triple for any W* satisfying Assumption 3.6. Since any operator Lp produced by the 
construction above for some choice of H can be produced using a different choice of distribu¬ 
tions but with H = I, one may conclude that each of these operators is part of a bispectral 
triple. (In other words, the obstruction to bispectrality that is visible when one seeks a bi¬ 
spectral triple for a given wave function ^ disappears if one instead focuses on the operator Lp 
and seeks a corresponding bispectral triple.) However, it seems plausible that some examples of 
bispectrality to be considered in the future will involve vacuum eigenfunctions that are neces¬ 
sarily non-commutative (unlike these examples in which the non-commutativity of the vacuum 
eigenfunctions can always be eliminated through a change of variables), and the observations 
and results above will prove useful in those contexts. 
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